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With much higher sensitivities due to coherence effects, it is often assumed that the first evidence 
for direct dark matter detection will come from experiments probing spin-independent interactions. 
We explore models that would be invisible in such experiments, but detectable via spin-dependent 
interactions. The existence of much larger (or even only) spin-dependent tree-level interactions is 
not sufficient, due to potential spin-independent subdominant or loop-induced interactions. We find 
that, in such a way, most models with detectable spin-dependent interactions would also generate 
detectable spin-independent interactions. Models in which a light pseudoscalar acts as the mediator 
seem to uniquely evade this conclusion. We present a particular viable dark matter model generating 
such an interaction. 



I. INTRODUCTION 

The sensitivity of dark matter (DM) direct detection 
experiments is undergoing rapid progress and is expected 
to continue in the next decade. There are a number 
of proposed experiments which will probe complemen- 
tary aspects of dark matter properties with much bet- 
ter sensitivities than the existing ones: DM mass, spin- 
independent (SI) and spin-dependent (SD) cross sections, 
the dependence of the cross sections on the target nuclei, 
directional information, etc. 

The focus, rightly, is often on the detection of spin- 
independent DM interactions, because, due to a coher- 
ence effect, the SI interaction cross section with heavy 
nuclei is enhanced by A 2 , the number of nucleons in a 
nucleus, and is, therefore, expected in many models to 
be the dominant interaction in DM detectors. 

There is a good chance that, in the not-too-distant 
future, direct detection experiments will be able to ex- 
tend their sensitivity to cover the full detectable param- 
eter space for SI cross sections, down to the 10 -48 cm 2 
level, below which atmospheric neutrinos constitute an 
irreducible background. 

Prior studies ^-Ql have considered the relationship be- 
tween SI and SD cross sections, concluding that the two 
are typically correlated when a viable dark matter candi- 
date is present. Most of the discussions have been in the 
context of the minimal supersymmetric standard model. 
(Similar statements have been made about DM candi- 
dates in universal extra dimensions jlj and little Higgs 
models ||, as well.) In general, the common wisdom is 
that SI experiments have a much better chance of first 
direct detection discovery. 

The generality of this conclusion cannot be addressed 
by merely considering operators; one must explore the 
underlying models which determine relationships be- 
tween operator coefficients. For example, the conclusions 
stated above ultimately stemmed from the assumption of 
DM with electroweak charges, which generically implies 
both mediators with at least weak-scale masses to justify 
null results thus far and couplings to the Higgs leading to 
SI signals. Once this condition is relaxed, the relation- 



ship between SI and SD cross sections becomes weaker, 
and models in which SD interactions are more easily ac- 
cessible, or even the only interaction accessible in direct 
detection experiments, become feasible. 

Here we point out that in order to impose the last 
condition, i.e., uniquely SD detection, the consideration 
of subleading effects is crucial. Since, due to coherence 
effects, SI experiments are more sensitive than SD ones 
(currently by 5 orders of magnitude), a loop-induced SI 
process might be only marginally more difficult or pos- 
sibly even as easy to detect than a tree-level SD one. 
Upon considering these additional operators, we find that 
models with light pscudoscalars are uniquely capable of 
generically evading such detection modes. 

Although several ingredients of our analysis appear in 
the literature || [|, the impact of light mediators on a 
general analysis of operators has not been heretofore dis- 
cussed, and the effect of loop corrections on DM scatter- 
ing has not been considered in this context. In Sec. 0, we 
review current bounds on SI and SD cross sections and 



the expected improvements. Sec. [II then constitutes the 



bulk of the paper. We discuss operators relevant for the 
detection of DM particles, including operators which be- 
come important in the case of light mediators. We then 
consider which models could generate exclusively SD in- 
teractions, and calculate the loop-induced interactions 
that would simultaneously be present. In Sec. [V we 



construct a viable model achieving our goals, in which 
the SI interaction is out of reach, but the SD interaction 
may be detected in future experiments. Sec. |y| concludes. 



II. PROSPECTS OF DIRECT DETECTION 

The best SI bounds come at present from 
XENONIO §, CDMS @, and XENONI00 @, with the 
highest sensitivity from XENONI00 near 3 x ICT 44 cm 2 
at 50 GeV. In general, optimal sensitivity is for DM 
masses of order the mass of the recoiling nucleus. At 
higher masses, the sensitivity decreases roughly as 
l/niDM- Within the coming years, XENON100, LUX, 
and SupcrCDMS can improve these bounds down to 
the 10 -45 cm 2 or possibly near the 10 -46 cm 2 level. 
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Ultimately multi-ton Xenon or Germanium experiments 
can achieve sensitivities to 10~ 47 cm 2 or maybe even 
10 -48 cm 2 , at which point atmospheric neutrinos 
form an irreducible background [^0| -^2|, and achieving 
sensitivity to lower SI DM-nucleon interactions seem 
unfeasible. 

For SD detection, the best current limit for DM-proton 
interaction is near 2 x 10 -38 cm 2 from SIMPLE 



slightly weaker bounds from COUPP 14 
and PICASSO ' 
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KIMS 



with 



|16|j , at similar optimal masses as above. 
For DM-ncutron cross sections the best bound is from 
XENON10 @ at 5 x 1CT 39 cm 2 at optimal sensitivity 
near 30 GeV. Within the next few years, COUPP @, 
PICASSO p), and XENON100 should improve these to 
few x 10~ 4 ^m 2 , for both protons and neutrons. These 
limits could then be extended to near 5 x 10~ 41 cm 2 with 
experiments such as DMTPC, or to 5 x 10~ 43 cm 2 for a 
500 kg extension of COUPP |o). 

Bounds on direct detection cross sections can also come 
indirectly from other experiments. One source is from 
DM annihilation signals from the Sun. The annihilation 
at equilibrium is proportional to the rate of DM cap- 
ture, which is driven by the same interactions as direct 
detection. In this case the SI terms are not so strongly 
enhanced over the SD ones, since this capture is mostly 
due to light nuclei, almost entirely hydrogen and helium. 
(Some enhancement does occur due to small amounts of 
Fe and O, but bounding the SD interaction neglecting the 
SI contribution is conservative.) Super-Kamiokande pi} ] 
and IceCube |2^] used this to place limits on SD proton 
interactions at around 10 -38 cm 2 , assuming annihilations 
primarily to bb. Above todm ~ 250 GeV, IceCube could 
even place a bound at 2 x 10 -40 cm 2 if the DM annihi- 
lated to W + W~ . However, these indirect bounds do not 
apply in the case of light mediators which will be dis- 
cussed below, since if the annihilations proceed through 
a light on-shell particle, decays to neither heavy quarks 
nor W bosons occur. 

Other bounds can be placed from constraints on op- 
erators from collider searches |23-25|. In cases where 
the mediator can be integrated out, these searches place 
bounds on interactions of very light dark matter better 
than those of direct detection, while remaining compet- 
itive with them for SD interactions of DM that can be 
directly produced at the Tevatron. The expected LHC 
reach is expected to also remain competitive with di- 
rect detection sensitivites of upcoming expierments [p4| . 
However, for mediators light enough to be produced on- 
shell, the bound deteriorates rapidly |25 , and is also not 



applicable for the class of models we discuss below. 



III. GENERAL CONSIDERATIONS 

A. Operator Analysis 

In order to survey possible models, we first identify 
all operators through which dark matter may interact 



Operator 


SI / SD 


Suppression 




= 4> 2 qq 


SI 






= (t> 2 qj 5 q 


SD 


q 2 


©3 


= ^d^q^q 


SI 




oi 




SD 


v 2 



TABLE I: Operators relevant for scalar dark matter detection. 
The suppression factor given is for the relevant cross section. 
Operators C9| and C| are only allowed for complex scalars. 



Operator 


SI / SD Suppression 


0{ = xx qq 
oi = xrfxqq 
oi = xxqW'q 
o{ = X7 5 x qi 5 q 


SI 

SI q 2 
SD q 2 
SD q 4 


oi = xn^xqip-q 
o{= x^r'xqj^q 

Oi = X7 m X97m7 5 9 
Ol = X7 M 7 5 Xg7M7 5 9 


SI 

SI v 2 
SD q 2 
SD v 2 or q 2 
SD 


oi = x^xq^^q 
o f 10 = xi^^xqo^q 


SD 

SI q 2 



TABLE II: Operators relevant for fermionic dark matter de- 
tection. Operators Og, 0( , Og, O{ only exist if the dark 
matter is Dirac. Notations as in Table \A. 



with detectors. In doing so, we will see which interac- 
tions give us the signals we are looking for, and which 
operators need to be suppressed by small coefficients or 
forbidden by symmetries. Similar operator analyses have 
been considered before in Refs. j| g, |2q| . We present it 
here as a guide to possible types of underlying structure. 

We assume that the mediator is heavy enough so that 
for the purposes of direct detection, describing the in- 
teraction of dark matter via a contact term is a reason- 
able approximation. Beyond this, we want to consider 
interactions with dark matter of arbitrary spin, without 
making additional assumptions, such as parity conserva- 
tion. At the structural level of the operators this en- 
compasses both clastic and inelastic scattering. Having 
two (or more) fields of different mass in the DM sector 
only leads to differences in kinematics and the presence of 
operators that are otherwise zero for Majorana fermions 
and real bosons for symmetry reasons (discussed below) . 

The smallest number of operators, as expected, are 
furnished by scalar dark matter candidates. These are 
listed in Table Q. Note that C| and C| are nonvanishing 
only if the dark matter candidate is complex. 

For fermionic dark matter, the operators are listed in 
Table |l| If the dark matter candidate is a Majorana 
fermion, the operators 0(, 0(, Og, O{ arc absent, 
as they are odd under charge conjugation. There are 
only two operators with tensor couplings. Since a^j 5 = 
i e v P a ( jp Cr /2, not all (pseudo)tensor-(pscudo)tensor com- 
binations are linearly independent. In addition, Oj has 
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Operator 


SI / SD 
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B^B,, qq 


SI 
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SD 
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01 = 


Bld»B»q lv q 


SI 
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Bld v B» g 7 „ 7 5 g 


SD 


v 2 
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B^d„B" q lv q 


SI 


v 2 q 2 
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B»d^B v q lvl 5 q 


SD 
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e^ pa B»d v BPq^q 


SI 
SD 


v 2 

<r 


oi = 


^ vpa B^d"B''q^q 


SD 





TABLE III: Operators relevant for vector dark matter detec- 
tion. Operators O3 and OI only exit for complex vectors 
fields. Notations as in Table fl. 



operators were suppressed or forbidden for some symme- 
try reason, the traditionally subleading operators could 
lead to contributions of the correct magnitude to be ac- 
cessible to current or future direct detection experiments. 
As pointed out in Ref. [^9) , these two statements may in 
fact be connected, since the spontaneous breaking of a 
symmetry forbidding the appearance of certain operators 
can provide for a natural explanation for the presence of 
light (pscudo) Nambu-Goldstone scalars. 

This opens up new possibilities. If SI operators with- 
out kinematic suppression factors are forbidden or highly 
suppressed for other reasons, the set of operators which 
may lead to a detectable SD signal becomes much larger. 

B. Renormalizable Models 



separate SD terms suppressed independently by v 2 and 
q 2 , while Og, commonly referred to as the anapole mo- 
ment coupling, has contributions to both SI and SD cross 
sections with different suppression factors. (Here, as else- 
where in the paper, v is the velocity of DM in the halo, 
approximately 10 -3 , while q is the momentum transfer 
in the interaction.) 



Finally, in Table III, we give the possible operators for 
vector dark matter candidates. Similar to the case of 
scalar dark matter, the operators O3 and 0\ arc only 
present if the vector is complex. 

There are a large number of operators that could medi- 
ate SD interactions. However, for our purposes, some of 
these may be ignored right away. For example, Og and 
O7 lead to both SD and SI interactions of comparable 
magnitudes. It may naively seem that all operators that 
come with kinematic suppression factors can be dismissed 
just as easily. After all, with DM in the galactic halo at 
such low velocities, the nonrelativistic limit is appropri- 
ate for detection, and traditionally such operators have 
indeed been neglected. Let us examine this assumption 
more carefully. 

Within the dominant WIMP paradigm, the mediator 
has typically been assumed to be at the weak scale, with 
direct detection occurring with 0(100 McV) momentum 
transfers and 0(100 kcV) recoil energies. In that case, 
the integrated-out mediator sets the scale of the oper- 
ators through a factor of In the nonrelativis- 
tic limit, terms like ?/>7 5 ?/; are suppressed by factors of 
\q\/2m,N or \q |/2mDM, while others, like V^^ 5 ?/;, have 
some components scale as v. Operators with any of these 
factor can typically be dismissed, because they are sup- 
pressed by O(10 3 ). This means that, even if present, such 
interactions can be ignored. For example, in the case of 
Majorana fermion dark matter, such as the neutralino 
in supersymmetric models, the only two operators that 
need to be considered are scalar - scalar and axial- vector - 
axial- vector jll|, |27], ; all others are highly suppressed. 

However, recent interest in explaining various possi- 
bly DM-relatcd anomalies have introduced models with 
O(GcV) mediator particles. In this case, if the leading 



If we wish to remain agnostic about the nature of the 
DM-nucleon interactions, we can say no more. However, 
if a further step is to be taken, it seems most conser- 
vative to assume that the DM comes from some theory 
with renormalizable interactions in which the operators 
leading to direct detection come from heavy states that 
have been integrated out. One can then ask what sort of 
renormalizable interactions could lead to the operators 
given above. Such a procedure was followed in Ref. [|J. 
Here we quote their results, along with the additional 
possibilities afforded by interactions yielding kinemati- 
cally suppressed operators. 

For scalar DM, the only option for generating solely 
SD operators seems to be a i-channel exchange of a light 
pseudoscalar, which yields 0|. While such an interaction 
breaks parity, given that parity is badly broken already 
in the standard model (SM), this is not a serious concern. 

For fermionic DM, several possibilities present them- 
selves. Once again a t-channcl light pseudoscalar ex- 
change produces solely SD interactions via O4. Addi- 
tionally, for Majorana fcrmions, the t-channel exchange 
of a vector with axial couplings, cither the SM Z or a new 
Z' , will generate only a single kinematically unsuppressed 
operator, Og . Other options are an s or u-channel cou- 
pling through either a scalar or vector, provided the cou- 
plings are chiral, in which case Og is generated again. If 
the couplings are not chiral, 0[ is produced as well. 

Finally, for vector DM, a light pseudoscalar in the t 
channel produces only O2, which breaks parity as in the 
scalar case. Alternatively an s or it-channel coupling 
through a fermion makes Og the leading operator, if the 
coupling is chiral while the vector boson is real. 



C. Loops and Subleading Interactions 

Suppose that one is presented with a model in which 
one of the above SD interactions is the only one present, 
or dominant over other by many orders of magnitude. 
Does that mean that only an experiment sensitive to SD 
interactions would see a signal? Not necessarily. 
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(a) Tree level 



(b) Loop processes 



(a) Tree level 



(b) Loop processes 



FIG. 1: The tree and loop level contributions to scattering of 
Majorana fermions through a Z boson. For all box diagrams, 
the crossed box diagram is included in calculations but not 
depicted. In the last diagram, a Higgs mediates the scattering 
through a Z loop. 



The bounds on SI cross sections are currently 5-7 or- 
ders of magnitude higher than the SD ones, and this looks 
to continue to be the case in the future. Therefore if any 
of the SD interactions discussed above induce subleading 
SI couplings, such an effect could potentially be visible 
in a SI experiment. There are two sources for such ef- 
fects. First, there are kinematically suppressed contri- 
butions of tree level scattering that were ignored above. 
These are easily estimated from Tables | |ni] given earlier. 
Second, the tree-level SD interactions can induce SI cou- 
plings at loop level. These are not as simple to estimate, 
and should be calculated to confirm their effect. 

Let us consider a Z (or Z 1 exchange) with a Majo- 
rana fermion, as in Fig. |la|. While the dominant contri- 
bution comes from 0(, also present is the anapole 
coupling. We see that this gives rise to a SI interaction 
suppressed by v 2 . Similarly, both the scalar exchange of 
Fig. ^a| and the equivalent diagram for vector exchange 
give an anapole coupling after using Fiertz identities. A 
fermion exchange of the same form in the case of vec- 
tor DM produces as well as Cg in the chiral limit, 
which again mediates a v 2 suppressed SI coupling. In all 
of these cases, there is a SI scattering cross section no 
more than O(10 6 ) smaller than the SD one, independent 
of any other field content of a model. This means that 
such interactions would be seen in SI experiments simul- 
taneously or in the next generation of experiements after 
they appear in SD ones. Only the pseudoscalar exchanges 
evade this, as they lead to no v 2 suppressed subleading 
contributions to DM-nucleon scattering at all. 

All the aforementioned interactions should also be 
computed at the one-loop level. While these will be sup- 
pressed by loop factors and extra couplings, they may 
also generate SI interactions. For large enough couplings, 
these loops might even give rise to interactions larger 
than the kinematically-suppressed ones discussed above, 
and so might be even more readily detectable. 

Without making any further assumptions about the 
underlying model, we can already identify diagrams 
which will produce SI interactions at loop-level. For SD 
interactions involving a i-channel exchange, at a mini- 
mum, exchanging two mediators in a box diagram will 
give rise to a SI interaction. For an s or u-channel pro- 
cesses, a SI loop level contribution can come from a loop 
with W or Z bosons exchanged between the quarks. 



FIG. 2: The tree and loop level contributions to scattering of 
Majorana fermions through a s-channel scalar. 



Consider the exchange of a Z with axial couplings to 
quarks. (We will discuss the case of a Z' shortly.) In 
that case, the quark level operator for tree-level scatter- 
ing (Fig. |la]) is 



9l 



rriQ Ti 1 - (J 5 - 5 

— T 3^—T XT 7 X 97^7 q , 
'w z m 7 



(1) 



where Q is the coupling of the DM to the Z. Then the 
DM-proton SD cross section generated is (see Apps. |a| 
and p| for details) 



SD 



(1.5 x 10" 39 cm 2 ) ( 

1 0.1 



(2) 



with the DM-neutron cross sections about 20% smaller. 
In this case, two one-loop processes lead to SI effective 
interactions: one with two Z exchanges, and a Higgs cou- 
pling through a Z loop to the DM (Fig. lb). We work in 
the limit m q <C mz *C uidm- (This limit is generally the 
one in which the DM has the correct relic abundance in 
models where the only coupling of the DM to the quarks 
is through electroweak bosons, while foregoing the last 
inequality only yields 0(1) changes, see Ref. [pOi.) The 
SI contribution to the effective coupling is then J3 0|, pl| 1 



9tQ 2 



4ir cos 4 9 



w mz 



(T§) 2 



2m| 



Am 2 h 



m q xx qq ■ 



(3) 



Taking a reference value of m/j = 120 GeV, this interac- 
tions will induce a SI cross section of 



< = (4 x 10- 47 



91 

0.1 



(4) 



Asking that the SD signal be just beyond current SD 
experimental bounds implies Q ~ 0.3, giving a SI cross 



1 In deriving this result, along with those following, we have set 
several quark operators, such as 



3m DM 

which all simplify to m q \X QQ on shell, but can have different 
nuclear matrix elements, to their on-shell value. In fact this 
seems to yield a conservative estimate, as out of the nuclear 
matrix elements known, the first one has the smallest value (for 
a detailed discussion of these issues see Ref. j33|). 
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v — V 

(a) Tree level (b) Loop processes 

FIG. 3: The tree and loop level contributions to scattering 
DM mediated by a light pseudoscalar. The dotted line can 
represent either a scalar, fermion, or vector boson. 

section of 4 x 10~ 47 cm 2 . This, while not detectable in 
experiments underway, is feasible with ones in prepara- 
tion. 

This result would make the v 2 -suppressed constribu- 
tion to SI scattering dominant. However, it is worth men- 
tioning that this cross section acts as a lower bound — it 
could be that the DM particle is part of larger represen- 
tation of 5(7(2), in which case additional loops involving 
Ws would also contribute. Generally, the size off the 
cross section grows as n 2 , with n the dimension of the 
representation |3l| , making it possible for the loop con- 
tribution to be dominant, and not merely competitive 
with the kinctically suppressed contribution, and even 
being large enough to be discovered simultaneously with 
the SD signal. 

If one wishes to consider models with a new Z' , then 
the existence of a Z' with Higgs coupling becomes model 
dependent. To talk about a lower bound, we can then 
ignore the contribution of the second term in the effec- 
tive coupling. The heavier mediator mass that such a 
model would entail would have to be offset with a larger 
coupling in order to be detectable. Thus, at loop level, 
one would generally expect the effective interaction to 
be of at least similar size, or possibly larger, due to 
the higher power of the coupling appearing in the loop- 
induced term. 

If one considers the possibility of a light Z 1 , which are 
not ruled out by collider constraints down to the GeV 
range for gauge couplings smaller than the SM by 10~ 2 , 
the situation discussed above would be reversed, and one 
would expect a smaller loop-induced contribution. How- 
ever, the SI contribution due to kinematically suppressed 
operators is insensitive to changes in the mediator mass, 
and would still be present. Constructing a model without 



such operators and without significant fine-tuning seems 
extremely difficult. It is difficult to say more in gener- 
ality, due to the large freedom in assigning masses and 
charges under a new gauge group. 

Now let us consider DM with chiral couplings to the 
SM via an s or u-channel. The most model-independent 
loop-level processes here come from box diagrams with 
the quarks exchanging a W or Z boson, an example of 
which is given in Fig. The contributions of the loops 
have completely different forms depending on whether 
the coupling of the DM is left- or right-handed. However, 
in all cases the loop-level processes only give rise to sup- 
pressed SD contributions. In addition to DM of the form 
in Fig. ||, this is also true for the cases of fermionic DM 
with a vector mediator and vector DM with a fermion 
mediator of similar topologies. In this case, we find that 
the most reliable lower bound on a SI cross section in 
this case comes from the v 2 suppressed contribution to 
the tree-level interaction discussed earlier. 

Finally, let us turn to the box diagrams induced in the 
cases of light pseudoscalar exchange, Fig. First, we 
consider the case of scalar DM. At tree level, the operator 
obtained after integrating out the pseudoscalar is 

-^o^Vqm^^^qi^q, (5) 
m 2 

where y q is the Yukawa coupling of the quark, so £ ab- 
sorbs both the coupling of the DM and mediator, and 
any scaling to Yukawas of the mediator-quark coupling. 
This leads to a a tree-level cross section of 

(See App. [b] for the definition of the cross section in cases 
of kinematically suppressed operators.) For a mediator 
with mass of a few GeV and £ = 0.01, this would be 
accessible to currently running searches. 

The calculation of the loop diagram in the same limits 
as the previous Z-mcdiated case does not give as com- 
pact of an answer, but can be expressed in closed form 
in terms of Passarino-Veltman scalar integrals |53| , com- 
puted with the use of FeynCalc Q as 



(47T) 



■evi 



C (m|,0,m|;m|,m^,0) - C {m\, m|, 0; m 2 , m\, m 2 ) 
f ml D (ml, ml, 0, 0, 0, m|; m 2 , m% ml, 0) ^d^^qj^q. 



(7) 



A numerical evaluation of the coefficients show the Cq and Dq functions with these parameters to scale as \n(m a /m^) 
and In 2 (m a /m^), respectively, beyond their overall 1/to? dependence. Using a fiducial value of m a /m,p = 0.01 gives 



i e-y 2 q 



(47T) 



, C s ^ 



(8) 
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where Cs ~ 80. Note that if the DM were real, this operator vanishes identically, and there is no loop induced 
coupling at one loop order at all. If present, the cross section induced is 



AN 



(4 x lO^ 4 cm 2 ) -2- 



£ \ 4 A00 GcV\ 4 



o.i J V ™t ) 



(9) 



undetectable for any any choice of parameters that would make the SD cross section detectable. 
The case of vector DM is very similar. For 



1 



(10) 



the tree- level cross section takes the same value as Eq. @. Meanwhile, the loop induced coupling is 
1 

(4^) 2 s y<? , 

+ Wo D (m 2 B ,m 2 B ,0, 0, 0, m|; to 2 , m 2 B ,m 2 a ,0) 



£ V q \ Co(m B ,0,m B ;m B ,m a ,0) - C o {m B ,m B ,0;m a ,m B ,m a ) 



4m 2 B 



B (m B ;m a ,m B ) - B o (m B ;0,m B ) 



(11) 



Bl d^B" 



which numerically evaluates to 



i evf 

(Ait) 2 m\ 



q CyB\ d"B u 



(12) 



with Cy ~ 80 very close to the scalar case, giving a loop induced SI cross section as in Eq. (^), and similarly giving 
no contribution if the DM were real. 

The case of fermionic DM is slightly different. This is because the tree-level operator responsible for scattering is 



(13) 



and, therefore, is parametrically suppressed by q , instead of the previous cases' q 2 . The tree level cross section then 
becomes 



ct& « (3 x 10- 43 cm 2 ) ( ^ 



nVi Gcv\ 4 



la J 



(14) 



We see that due to the greater momentum suppression, we require a lighter mediator mass and cannot afford the 
coupling of the DM to be as small as in the bosonic case above. In this case a cross section detectable in current 
experiments would require, for example, a mediator with m a = 100 MeV and £ = 0.1. 

Meanwhile, the effective coupling from computing the loop diagram in the same limits as the other cases is 



i evi 

(47I - ) 2 to 2 



1 777 

i X r< I 1 2n 2 2 2\ 2 r~i in 2 2 n 2 2\ 

c o(m x , m x , 0; m a ,m % , m a ) - m x C o (0, m x , m x ; 0, m a , m x ) 



XTXTfpQ 



1 + B a (m 2 x ;0, m 2 ) - B (0;m 2 a ,m 2 a ) + Am 2 , C (m x , 0, to 2 ,; m 2 , to 2 , 0) - m 2 C (to 2 , to 2 , 0; to 2 , to 2 , to 2 ) 



+ 3to 2 to 2 D (m x , m 2 , 0, 0, 0, to 2 ; to 2 , to 2 , to 2 , 0) 



XX 11 



(15) 



which numerically yields 

i evi 



(47r) 2 to 2 



Cf 1 Xl^X Qltil + Cf 2 XX QQ 



(16) 



with ~ 4.8 and Cf 2 ~ 170. The magnitudes of 
these coefficients can be understood as arising from the 



ln(m /m^) and ln 2 (TO a /TO^) behavior of Co and D men- 
tioned above. The loop-level cross section is then 



<j$? « (3 x 10- 56 



£ \ 4 /100 GeV x 4 



(17) 



0.1 / \ TO X 

We will confirm below in the explicit model of Sec. IV 
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that the loop induced coupling is indeed tiny, but it is 
simple to see here why this is generically so. 

Unlike in the massive mediator cases, there are two 
mass scales in the dark sector, that of the DM itself and 
that of the mediator. At tree level, the lighter media- 
tor mass is the one that appears in the denominator of 
the operator. However, at loop level, the value of the 
loop integral is parametrically controlled by the mass 
of the DM, the heaviest particle in the loop. Addition- 
ally, a pseudoscalar which is the Nambu-Goldstone bo- 
son of a broken symmetry would be expected to couple to 
quarks proportional to the masses of the quarks. Thus, 
at loop level, the effective operator would be expected 
to be suppressed by extra factors of quark Yukawa cou- 
plings. Together, both effects combine to make the loop- 
level coupling to be as many as 20 orders smaller than 
the tree level one, with higher order corrections to the 
non-relativistic scattering approximation coming at sim- 
ilar orders as q 2 v 4 , so that the SI induced interaction is 
expected to be completely negligible. 

IV. THE AXION PORTAL 

We have just seen that without tuning of couplings, 
models with light pseudoscalar mediators provide the 
unique method of avoiding any SI signal, while still pro- 
ducing a SD direct detection signature. Now we turn to 
the question of whether a viable model producing DM 
with the correct abundance can have these features. 

Coupling a light pseudoscalar to quarks is most effi- 
ciently achieved by adding a scalar field which sponta- 
neously breaks a global symmetry, and which, by mix- 
ing with the Higgs, gets a coupling to the SM. Allowing 
this scalar to have a new global charge, while adding 
new fermions charged under the same symmetry, ensures 
that the new scalar field is the only method for the new 
fermions to interact with the SM. 

As a simple realization of such a mechanism, where 
the dominant interaction is O a via a pseudoscalar inter- 
action, we introduce, following Ref. |3f|, a scalar field 
charged under a new global U(l)x charge that is spon- 
taneously broken to 

s =(^Mik)- (18) 

This scalar field is coupled to a new fermion, which is 
vector-like under the SM, through C = — £5xX c + h.c., so 
that after the scalar field acquires a vacuum expectation 
value, the fermion receives a mass of m x = £/ a , allowing 
it to act as dark matter, with stability ensured by the 
remnant of U(l)x after breaking. 

In order for the pseudoscalar to interact with the SM, 
some known particles must also carry charges under the 
new (7(1) x- In a two Higgs doublet model, this can be 
accomplished by adding a term of the form 

C = XS n H u H d + h.c, (19) 



by assigning the appropriate charges to the Higgses and 
SM fermions, and promoting the U{l)x to a Peccei- 
Quinn (PQ) symmetry. For n — 2, this cou plin g is of the 
same form as in the case of the DFSZ axion |36| , |37| , while 
the n = 1 case functions like that of the PQ-symmetric 
limit of the next-to-minimal supcrsymmetric standard 
model |3^] We now have a dark matter candidate cou- 
pling to the SM though a massive scalar and an axionlikc 
Nambu-Goldstonc boson. The Nambu-Goldstonc boson 
is assumed to get a small mass through an unspecified 
mechanism. Anticipating making the scalar heavy, by 
virtue of 

m 2 / m 2 \ 

a choice of, say, m s = f a = 1 TeV and m x = 1.1 TeV 
(corresponding to £ = 1.1) yields a cross section of 
3 x 10 ~ 26 cm 3 /s and so generates the correct order of 
magnitude for the relic abundance |p5| . 

For direct detection, two channels present themselves. 
The scalar gives a SI cross section through the operator 
0{, due to mixing of the scalar with the two CP-even 
Higgses, while the light axionlike state yields a SD inter- 
action, C{, by a similar mixing with the CP-odd Higgs. 
For our purposes, we need check whether this tree-level 
SI cross section can be small enough to be completely 
negligible. 

The mixing of the scalar with the two CP-even Higgs 
has a lot of arbitrariness to it due to the 1 1 constants in 
the most general [/(l)pQ-preserving two-Higgs-doublet 
and one-singlet potential. However, we can say that 
barring accidental cancellations, this mixing will be e = 
0(y ew / f a ), so that we may write the tree-level SI cross 
section as 

cr*f « (2 x 10- 42 cm 2 K 2 £ 2 (^^1^ , (21) 

(See App. |A] for a caveat on the values of the nuclear 
matrix elements in this calculation.) In the model con- 
sidered in Ref. m s needed to be light, O(10 GeV), in 
order to provide a mechanism for Sommerfeld enhance- 
ment to explain astrophysical anomalies. In that case, 
the direct detection cross section was in tension with the 
SI bound and could only be slightly beyond current lim- 
its, at a few x 10~ 43 cm 2 . However, if we impose no such 
condition, m s could be larger. If it is at the electroweak 
scale, then the cross section is at most a few x 10~ 45 cm 2 , 
smaller than the sensitivity of the next generation of di- 
rect detection experiments. If m s ~ 0(1 TeV), a rea- 
sonable choice given the scale of f a in this setup, then 
the cross section becomes undetectably small, below the 
irreducible atmospheric neutrino limit. 

Let us next consider the pseudoscalar channel. With 
the interaction kinetically suppressed by the momentum 
transfer as q 4 , wc cannot merely compute the cross sec- 
tion in the limit of q 2 — > as we did in the scalar exchange 
case. Instead, we must define a cross section at a fixed 
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momentum transfer. (See App. [b] for a more thorough 
discussion.) We choose to do so at q 2 ei = (100 MeV) 2 . 
Because the signal is different from that of unsuppressed 
interactions relative to the expected recoil energies, the 
sensitivities of experiments are modified. This was stud- 
ied in Ref. |2!| , with the result that at the same reference 
momentum transfer, optimal sensitivities of SD experi- 
ments to pscudoscalars remained at the same order of 
magnitude as in the unsuppressed case, but with 1 / muM 
scaling of the limits. 

With this definition, we can compute the SD cross sec- 
tion for q 2 = q 2 cf as 



7 SD 



(2 x 10~ J ' cm^sirT 



4m 2 



1 GeV 



m a 



(22) 



where tan0 = n sin 2/3 [v cw /(2f a )] is the mixing of the 
s with the Higgses |35|. From this we see that given 
a DM mass m x = 1.1 TeV, a pseudoscalar with a mass 
m a w 300 MeV generates a cross section of 3 x 10 -40 cm 2 , 
within the range of the next generation of direct SD de- 
tection experiments. In fact, in a two Higgs doublet 
model like this, the nuclear matrix element also has a 
dependence on /3 as up-type quarks couple with a coeffi- 
cient proportional to cot /3, while down type ones couple 
proportional to tan/3. We have evaluated the matrix el- 
ements for the above cross section at tan/3 = 1. At large 
values of tan j3 the cross section can rise by almost 2 or- 
ders of magnitude. 

Given the tiny size of the tree-level SI cross section, and 
in keeping with the discussion of the previous section, we 
should confirm that the loop-induced couplings fail to 
produce a detectable SI cro ss sec tion. The calculation 
mostly mirrors that of Sec. Ill C . The only substantial 
difference is the aforementioned different coupling to up 
and down type quarks. As before, we evaluate the nuclear 
matrix elements at tan /3 = 1, but this time, varying tan /3 
cannot only modify the cross section by a factor of O(l) 
as the suppression of sin 2 9 at high tan /3 is too strong, so 
we find 



'si 



(3 x 10" 5b cm 2 ) 



£ sin6»\ 4 /100 GeW 4 
0.1 



(23) 



with no additional implicit tan/3 dependence. 



V. CONCLUSIONS 

As the sensitivity of both SI and SD direct DM de- 
tection experiments increases, it is worth asking to what 
extent the discovery potential of the two methods is com- 
plementary. In this work we have pointed out that when 
one considers the full range of possible mediators, instead 
of being confined to new weak-scale particles, the range 
of possible viable interactions generating SD cross sec- 
tions increases. At the same time, when one searches for 
interactions for which SD experiments are complimen- 
tary for discovery — ones which could not be seen in any 



SI experiments without the need for accidental cancel- 
lations or other tuning — it becomes necessary to take 
into account sublcading contributions to scattering, such 
as suppressed operators and loop processes. The outcome 
is that the traditional models considered also generically 
produce SI interactions whose suppression is counterbal- 
anced by the greater sensitivity of SI experiments. The 
list of viable candidates whose interaction with the SM 
can be described by tree-level mediators integrated out 
in a renormalizable model is then reduced to merely ones 
mediated by light pseudoscalars. 

We have presented a realistic model of such interac- 
tions that generates the right DM abundance with a 
fermionic DM candidate without having other interac- 
tions generating detectable SI interactions. 

Similar scenarios can also be considered with a scalar 
or vector dark matter candidate. Just as in the case 
of fermionic DM, O s { v gives the leading interaction in 
the nonrelativistic limit, while O^" is kinematically sup- 
pressed. The necessary couplings between the pseu- 
doscalar and the scalars or vectors cannot be generated in 
as simple a manner as those used above, so more model 
building will be required. However, the suppression is 
only by q 2 , so the mass differences between the scalar 
and pseudoscalar do not have to be quite as large, and 
the couplings themselves can be smaller, so that the pa- 
rameter space of couplings and the pseudoscalar mass are 
not as tightly limited by experiment, potentially making 
the exercise worthwhile. 
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Appendix A: Nuclear Matrix Elements 

Here we summarize how to compute the dark matter- 
nucleon interaction cross sections from quark-level inter- 
actions. Much of this has been discussed in the DM lit- 
erature, with the exception of the pseudoscalar matrix 
element, as it only plays a role in momentum suppressed 
cross sections. 

For a vector coupling, nuclear matrix elements are 
straightforward to compute, since a vector coupling to 
quarks is a conserved current, so the coupling to a nu- 
cleoli is obtained from the sum of the currents of the 
valence quarks. 

In the case of a scalar coupling to quarks, we are in- 
terested in the effective nucleon coupling induced by a 
quark level coupling: 



a q m q qq 



f N m N NN. 



(Al) 
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We define the nuclear matrix elements conventionally by 
(N\m q gg\N}=m N f^\ (A2) 

On including the coupling to gluons induced by integrat- 
ing out heavy quark loops, /jy is given by 



(AT) 



JTq 

q—u : d,s 



a ^ ( A3 ) 

q— c.b,t 



where ffQ 



(AT) 



1 En=ti,ii,s frq 

Unlike the u and a matrix elements, which can be ex- 
tracted from 7riV scattering, the uncertainty associated 

with the strange quark matrix element f^J is higher, 
which introduces a substantial uncertainty in the SI cou- 
pling to nuclcons. Most studies use numerical values 
frpP 3> fxud based on older calculations. A represen- 
tative set of values is that used by the DarkSUSY pack- 
age p0|, wherein, 



f- 



(p) 



Tu 



0.023. 



4 r : } =0.019, 



Ap) 

JTd 
An) 
JTd 



0.034, /<£ 



0.041 



,(n) 
JTs 



0.14, 
= 0.14. 



(A4) 



These are the values used for the numerical estimates 
given above, and in most of the literature. However, re- 
cent lattice QCD results give substantially smaller values, 
fW = 0.013 ± 0.020 @ (see also || ||), and so the 
SI cross section from scalar exchange (if it couples pro- 
portionally to mass) may be smaller by a factor of 2 - 5 
than numerical results quoted by many calculations. 

For SD interaction we need to consider the nuclear 
matrix elements induced by the quark level axial- vector 
and pscudoscalar couplings, 



a N NsWn, 



(A5) 



and 



c q m q qij 5 q — > g^ tun Ni'j N. (A6) 
For the axial- vector current, defining 

(N\q 7 ^ 5 q\N} = s^Aq( N \ (A7) 



where is the spin of the nucleon, we have 



q—u,d,s 



(N) 



(A8) 



The matrix elements coming from polarized deep inelas- 
tic scattering carry much smaller uncertainties than for 
the scalar SI interaction above. For our numerical results, 
we use again the DarkSUSY values, 

Aii( p » = Ad^ = 0.77 , 

Ad {p) = Au (n) = -0.40 , 

As^ = As (n) = -0.12. (A9) 



More recent determinations favor slightly different val- 
ues, and the PDG quotes As*") = -0.09, Ad^ = 0.84, 
Au (n) = -0.43, with a 0.02 uncertainty for each Q; the 
effect on our numerical results is neglig ible . 

For the pseudoscalar current in Eq. ( |A6| ) the nucleon- 
level coupling is determined by the same axial-vector 
matrix elements above. The relationship is estab- 
lished through generalized Goldbcrger-Treiman rela- 
tions. While not normally considered in dark matter 
detection, it has been well-studied in the axion litera- 
ture |45|, Q . Taking divergences of the axial currents and 
using the equations of motion for the quarks yields p^] 

9n = (c u - c q rj) Au (Ar) 

+ (c d - c qV z) AdW + (c s - c qV w) AsW , (A10) 

where i] = (1 + z + w) , z = ra u jra^ and w = m u /m Sl 
while c q is the mean of the quark coupling coefficients. 
Due to uncertainties in the value of z, the value of g^ 
can vary by as much as a factor of 2. 



Appendix B: Cross Sections 

In this Appendix, we provide a summary of cross sec- 
tions for DM-nucleon interactions relevant for calculating 
the various cross sections discussed above in the non- 
relativistic limit. 

Wc first consider the unsupprcsscd opcrors in the limit 
of zero momentum transfer. SI cross sections can come 
from either scalar or vector quark couplings. Effective 
DM-nuclcon scalar interactions for fcrmions of the form 



f N XxNN, 



(Bl) 



which are derived from the quark-level couplings using 
nuclear matrix elements, as explained in App. |A| lead to 
a DM-nucleus cross section 



■jP[Zf p + (A-Z)f n 



(B2) 



for Majorana DM fcrmions. (For Dirac fermions, all re- 
sults for Majorana fcrmions arc divided by 4.) Here fx is 
the reduced mass of the DM-nucleus system. The per- 
nucleon cross section, which is usually quoted for com- 
parisons, is 



a = pi 2 ^[Zf p + (A-Z)f n 



(B3) 



where /i is the reduced mass of the DM-nuclcon system. 

For scalar or vector dark matter, the relevant operators 
are (we include the DM mass to give all operators the 
same dimension) 



f N m4, # NN or f N m B B^B^NN, 
the nucleon cross section for either operator is 



1 



1 



V 2 ^[Zf P + (A-Z)f n 



(B4) 



(B5) 
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Vector interactions for fermions only exist in the case 
of Dirac DM: 



The only other unsuppresscd SD interaction is for vector 
DM and comes from 



(B6) 



where, b p 



2&u + bd and b n = b u + 2bd, due to vector 
current conservation, as discussed above in App. |A]. Then 



1 2 1 

7T A Z 



^■[Zb p + (A-Z)b n y. (B7) 

For the operators 

bN^d^N^N or b N B^d^B" N^N, (B8) 

which only exist for complex scalars or vectors, the cross 
section is 



1 



1 



a= -M -To [Zb p + (A-Z)b r , 

7T A Z 1 



(B9) 



Unsuppressed SD interactions come solely from the 
quarks' axial currents. In the case of 



a N Xl^xNj^N, 



the DM-nucleus cross section is 

o = — fi 2 a 2 N J N (J N + 1), 

7T 



and for a nucleon 



2 2 

a = — fi a N . 

7T 



(BIO) 



(Bll) 



(B12) 



a N e^PB^Bv N lpl 5 N . 



(B13) 



Here, the DM-nucleon cross section is 



2 2 

a = — /i a N . 

7T 



(B14) 



All of the above cross sections are quoted in the q 2 — > 
limit. In this limit, interactions mediated by light pseu- 
doscalars are all zero, so we need another way of ex- 
pressing such cross sections. To do so, we will use the 
fact that while in the nonrelativistic limit tpip ~ 2m, 
■07 s -0 ~ q l ^a 1 ^, so that using the results above we can 
write (since q 2 ss \q\ 2 in the nonrelativistic limit). 



4m| 4m^ 



XI^XN^N . (B15) 



We then compute the cross section as above, and quote 
a result at a reference value of q 2 . We have chosen 
q 2 = (100 McV) 2 since with q 2 = 2m N E R , where E R 
is the recoil energy of the nucleus, this is a typical value 
for most SD detectors. Other momentum-suppressed op- 
erators can be handled the same way. 
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